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In this paper, we study theoretically and experimentally the friction between a rough parabolic or conical 
profile and a flat elastomer beyond the validity region of Amontons' law. The roughness is assumed to be 
randomly self-affine with a Hurst exponent H in the range from 0 to 1. We first consider a simple Kelvin 
body and then generalize the results to media with arbitrary linear rheology. The resulting frictional force as 
a function of velocity shows the same qualitative behavior as in the case of planar surfaces: it increases 
monotonically before reaching a plateau. However, the dependencies on normal force, sliding velocity, shear 
modulus, viscosity, rms roughness, rms surface gradient and the Hurst exponent are different for different 
macroscopic shapes. We suggest analytical relations describing the coefficient of friction in a wide range of 
loading conditions and suggest a master curve procedure for the dependence on the normal force. 
Experimental investigation of friction between a steel ball and a polyurethane rubber for different velocities 
and normal forces confirms the proposed master curve procedure. 

Friction is a phenomenon that people have been interested in for thousands of years but its physical reasons 
are not clarified completely yet. Not only is it still not possible to predict the frictional force theoretically, 
there are also no reliable empirical laws of friction which would satisfy the needs of modern technology. In 
practice, the simplest Amontons' law of dry friction 1 is usually used, stating that the force of friction, F, is 
proportional to the normal force, F N : F = /.iF N , the proportionality coefficient ji being called coefficient of friction. 
According to Amontons, the coefficient of friction does not depend on the normal force and the contact area. 
Amontons did not differentiate between the static and the sliding coefficients of friction, nor even between 
different materials (he states that the ratio of the frictional force to the normal force is "roughly" one third of 
the normal force, independently of the contacting materials as long as they are not lubricated 1, pp " 20S_209 ) . However, 
already Coulomb knew that the coefficient of friction, even between the same material pairing, can change by a 
factor of about four depending on the contact size 2, p 43 and on the normal force 2, p 66 and B . As a matter of fact, there 
are no obvious reasons for the validity of Amontons' law. On the contrary, much effort has been made in the 
1940th-60th years to understand why Amontons' law is approximately valid 3,4 . More recently, the violations of 
Amontons' "law" became a subject of vivid interest among researchers and focus has shifted towards a detailed 
understanding of factors responsible for these violations as well as the construction of generalized laws of friction. 
Several recent studies focused on violations of Amontons' law for static friction due to the dynamics of onset of 
sliding, which have been studied both experimentally 5 7 and theoretically in three-dimensional 8 and one-dimen- 
sional models 9,10 . These works showed that it is possible to understand and to exactly control the "violations" of 
the Amontons' law. The detailed dependence of the sliding coefficient of friction on normal force was not studied 
yet. In the present paper, we investigate the force of sliding friction beyond Amonton's law for elastomers with 
linear elastic rheology, and formulate rules for constructing generalized laws of friction for this class of materials. 
Theoretical findings are supported by experimental data. 

We use the standard assumption that elastomer friction is caused by energy dissipation in the volume of the 
elastomer due to local deformations 11,12 . For illustrations and basic understanding we implement this mechanism 
in the frame of a one-dimensional model based on the method of dimensionality reduction (MDR) 13 . As this 
method has recently been subject of a very controversial discussion, we would like to briefly describe the 
arguments for and against its validity. For small normal or tangential movements of three-dimensional bodies 
in contact, the total energy dissipation can be calculated either by integrating the dissipation rate density over the 
whole volume or equivalently determined over total force and displacement. Either method, providing the correct 
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relation between macroscopic forces and displacements is therefore 
suited for simulation of elastomer friction. It was shown that the 
force-displacement relation of elastic bodies can be described by a 
contact with a properly defined one-dimensional foundation, pro- 
vided the indenter is an (arbitrary) body of revolution 13 . As the 
macroscopic contact configuration is one of the essential properties 
in our treatment, we conclude, that the MDR is appropriate to simu- 
late this "macroscopic part" of the contact problem. Furthermore, 
the applicability of the method of dimensionality reduction (MDR) 
(again with respect to the force-displacement relations) was later 
extended and verified for contacts of randomly rough, self-affine 
surfaces with elastic 14 and viscous 15 media and for tangential con- 
tacts 16 . In the paper by Lyashenko et. al. 17 , the validity of the MDR 
was questioned. The authors of this paper showed that the MDR does 
not allow the simulation of the real contact area of random surfaces 
(apart from the limit of very small forces). This conclusion is correct 
and is confirmed by Popov in the review paper 18 and the mono- 
graph 16 devoted to the MDR. The force-displacement relations, on 
the contrary, are correctly described in the frame of the MDR. This 
was illustrated in 16, cha P ter 10 by comparison with direct three-dimen- 
sional simulations of randomly rough surfaces having the spectral 
power density C(q) ~ q~ 2H ~ 2 , where q is the wave vector, for values 
of parameter H in the range — 1 < H < 3 as well as for contacts of 
rough spheres 19 . This shows, that the MDR can be a useful modeling 
tool even for simulation of fractal rough surfaces as long as we stay 
within the standard Grosch' paradigm of the rheological nature of 
elastomer friction 12 . 

Following the above arguments, we use the MDR as a basis for the 
simulations of elastomer friction in this paper. As our aim is achiev- 
ing a basic understanding of frictional behavior between a rigid 
body and an elastomer, we consider the following simple model: 
(a) the elastomer is modeled as a Kelvin body, which is completely 
characterized by its static shear modulus G and viscosity }], (b) the 
non-disturbed surface of the elastomer is planar and frictionless, (c) 
the shape of the rigid counter body is assumed to be a superposition 



of a simple macroscopic profile and random self-affine fractal 
roughness without long wave cut-off (Figure 1), (d) no adhesion 
or capillarity effects are taken into account, and (e) we consider a 
one-dimensional model. In the subsequent sections, the restrictions 
are discussed and the model is generalized for arbitrary linear 
rheology. 

Results 

Generalized law of friction for rough indenters with power-law 
profiles. Let us consider a rigid indenter having the form 

z=g(x)=g 0 (x)+h(x) (1) 

consisting of the macroscopic power-shaped profile 

go(x) = c„\x\" (2) 

and a superimposed roughness h(x), as shown in Figure 1. 

Coordinates x and z are measured from the minimum of the 
macroscopic form, so that g 0 (0) = o. The ensemble average of the 
rough profile is assumed to be zero: (h(x)) = 0. The roughness was 
assumed to be a self-affine fractal having the power spectral density 
Cm oc q~ 2H _1 , where q is the wave vector and H, the Hurst exponent. 
This one-dimensional power density corresponds to the two-dimen- 
sional power density of the form C 2 d oc q~ 2H ~ 2 1S . The spectral 
density was defined in the interval from q min = 2n/L, where L is 
some reference length, to the upper cut-off wave vector q max = n/ 
Ax. The spacing Ax determines the upper cut-off wave vector and is 
an essential physical parameter of the model. Surface topography was 
characterized by the rms roughness h 0 , which is dominated by the 
long wavelength components of the power spectrum and the rms 
gradient of the surface Vz, dominated by the short wavelength part of 
the spectrum. Throughout this paper, we will assume that the 
indentation depth of the indenter, d, is much larger than the rms 
value of the roughness, ho<$id. This means that the large-scale con- 
figuration of the contact is primarily determined by the macroscopic 
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Figure 2 (a) A "large scale picture" of a contact of an elastomer and a rigid conical indenter which is moving with velocity v. (b) Rheological model for a 
viscoelastic medium. 



form of the indenter and does not depend on the roughness 
(Figure 2a). 

According to the MDR 18 , the Kelvin body can be modeled as a 
series of parallel springs with stiffness Ak z and dash pots with damp- 
ing constant Ay (Figure 2b), where 



Ak z = AG Ax, Ay — 4r\Ax. 



(3) 



The rigid indenter is pressed into a viscoelastic foundation to the 
depth d and is moved in tangential direction with velocity v 
(Figure 1), so that at time t it is described by the equation 

z=g(x + vt)-d=g(x)-d. (4) 

For convenience, we introduced the coordinate x = x + vt in the 
coordinate system moving together with the rigid indenter. 

The normal force in each particular element of the viscoelastic 
foundation is given by 

/= -4Ax[Gu(x)+riu(x,t)] (5) 

where u is the vertical displacement of the element of the viscoelastic 
foundation. For elements in contact with the rigid surface, this means 
that 

f = 4Ax[G(d-g(x))-rivg'(x)]. (6) 
The normal and the tangential force are determined through equations 



F N =4 



[G(d-g&))-TivgfQc)]d3c, (7) 



g'(x)[G(d—g(x))—rivg'(x)]dx. (8) 



We first consider the force of friction at very low velocities. The contact 
configuration is then approximately equal to the static contact. The 
uppermost left and uppermost right points — <X\ and a 2 of the contact 
(see Figure 2a) are then both determined by the condition g{— — d 
~ g( a i) ~ d = 0. Because of the relation ^—aj) = gia 2 ), the integrals 

«2 " 2 



t]vg'{x)Ax and 



Therefore 



F N = 4 



G{d-g(x))dx=4 



G(d-g 0 (x))dx, 



(9) 



F x =4r)v [g'ixfdx^nv [g 0 (x) 2 + h(xf]dx (10) 



We assume that the gradient of the macroscopic shape of the indenter 
is much smaller than that of the roughness, (g 0 (x) 2 ) <C ( h (x) 2 ) , so that 



F x ^4i]v 



h (x) 2 dx«4;jvVz 2 L C( 



where Vz is the rms value of the surface gradient and L co 
the contact length. For the coefficient of friction, we get 

4L cont v)i 2 4GL cont vz 2 
- Vz = Vz 



Fn 



F N 



(11) 



+ a 2 



(12) 



where t = }]IG is the relaxation time. This equation shows, that both 
the macroscopic shape of the indenter and the microscopic properties 
of surface topography determine the coefficient of friction: the contact 
length is primarily determined by the macroscopic properties (shape of 
the body and the normal force) while the rms gradient is primarily 
determined by the roughness at the smallest scale. 

Consider the opposite case of high sliding velocities. The detach- 
ment of the elastomer from the indenter occurs when the normal 
force (which is the sum of elastic and viscous force, Eq. (6)) vanishes. 
If the rms value of the elastic force and the viscous force become of 
the same order of magnitude, the detachment will occur in almost all 
points with negative surface gradient, thus a one-sided detachment of 
the elastomer from the indenter will take place. Characteristic rms 
values of the three terms in Eq. (6) are proportional to Gd, Gh, and 
f/vVz. If the indentation depth is much larger than the roughness of 
the profile, d^>h, then the condition for one-sided detachment of the 
elastomer from the indenter reads Gd ~ qvVz or 



d/(TvVz)«l. 



(13) 



In that case, the friction coefficient achieves an approximately con- 
stant value 20 of 



= v/2Vz 



(14) 



Although this result was derived in ref. 20, we would like to comment 
here shortly on it. In the plateau region, the elastomer behaves prac- 
tically as a viscous fluid: the elasticity does not play any role and all 
contacts are "one-sided." The normal and tangential forces reduce to 



g'{x)G{d-g(x))dx in (7) and (8) vanish. F x = 4 



r\v{z'{x)) 2 dx, F N =4 



(real cont) 



(real cont) 



// v | z f (x) I dx. For the 



normalized coefficient of friction we get 



_ falcon,) Wgd* _ v (/(real co„„ 

^~ ^)|d7^ VZcont 



Ji 



(real cont) 



1/2 



/(real cont) W<.*)\** 



(15) 



For an exponential probability distribution function of the gradient 
of the surface, the ratio of the integrals in (15) is equal to y/2, in 
accordance with (14), and it depends only weakly on the form of the 
distribution function. 
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Figure 3 | Dependence of Ji on £ at low velocities. Different symbols 
correspond to different sets of parameters H, v, F N , G, 8, t and Vz. All data 
collapse to one master curve described by equation (17). 

For the macroscopic power law shape (2), the indentation depth 
and contact radius are given by 16 



8Gc„n 



F N («+i) c y»v +i 



8Gn 



(16) 



Substituting the contact length L cont = 2a into equation (12), we 
obtain the coefficient of friction at low velocities: 



where we introduced dimensionless variables 
vtVz 



v/2Vz 



V2 



n + l 

c„n 



17) 



(18) 



Numerical simulations presented in Figure 3 for the case of a rough 
cone show that all data in the coordinates collapse to one 

master curve with a slope equal to one. The validity of equation 
(17) was numerically confirmed for the following ranges of para- 
meters. The reference length of the system was L = 0.01 m and 
the number of elements N — L/Ax was typically 5000. 11 values of 
the Hurst exponent ranging from 0 to 1 were studied. All values 
shown below were obtained by averaging over 200 realizations of 
the rough surface for each set of parameters. Parameter studies have 
been carried out for 20 different normal forces F N ranging from 10~' 
to 10 1 N, 20 values of the G modulus from 10 5 to 10 7 Pa, 20 values of 
rms roughness h 0 from 10~ 6 to 10~ 4 m, and 20 values of the spacing 
Ax from 10~ 7 to 10~ 5 m, 20 values of angles 0 ranging from 5 to 75°, 
and 20 relaxation times z ranging from 10~ 4 to 10~ 2 s, while in each 
simulation series only one parameter was varied. 

It is easily seen that we can get both limits (12) and ( 1 4) by writing 

i=( 1+r '^ (19) 



where a is a dimensionless fitting parameter. Numerical simulations 
(Figure 4) show that this dependence is valid for all parameter sets 
used in our simulations, while the best fit is achieved with a = 1.5. 
Interestingly, parameter a seems not to depend on the macroscopic 
shape of the indenter. 

Friction law in a case of a general linear rheology and the "force 
master curves". Let us now consider friction of elastomers with a 
more realistic rheology, which is characterized by the frequency 
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Figure 4 | Dependence of jx on if for the following set of parameters: F N = 
10 N, G= 10 6 Pa, h 0 = 5 • 1(T 5 m, r = 10" 3 s, H= 0.4 and 0= 10° 
(conical indenter) or R = 10~ 2 m (parabolic indenter). The solid line 
corresponds to the analytical approximation (19) with a = 1.5. 

dependent complex shear modulus G(co) = G'(co) + iG"{(D), 
where G' is the storage modulus, and G" the loss modulus 21 . At 
low frequencies, the shear modulus tends towards its static value 
G 0 - For simplicity, we will assume that the macroscopic contact 
mechanics of the indenter is completely governed by the static 
shear modulus G 0 , which is correct for sufficiently small sliding 
velocities. On the other hand, the frictional force is almost 
completely determined by the smallest wavelength components in 
the spectrum of roughness and thus by high frequency rheology. The 
frictional force at low velocities can be therefore estimated by using 
equation (11) and substituting — » G"(<x> m ax)/ (B max> where 



. 4G"(C0 max )vVz 2 L cont 4G"(«?maxV)V Z 2 I c , 



(20) 



The contact length is completely determined by the macroscopic 
contact mechanics of the indenter, equation (16), where we 
substitute the constant static shear modulus G 0 . For the friction 
coefficient, we therefore get 



4G"(g max v) {F N (n + l 



Vz 2 



(21) 



At high frequencies, the plateau value of will be achieved. An 

interpolation between (21) and this value is provided by 

-1/7 



v/2Vz 



H 1 + 



8G 0 c„n 
2x/2G"(q ma xV)Vz \F N {n + l 



(22) 



with a. ~ 1.5. This equation shows that the coefficient of friction for 
rigid bodies having macroscopic power-law shape has the general 
form 



H~Q(F N -p(v)), 



(23) 



where p( v) is a function of velocity, which depends on the rheological 
properties of the elastomer. Since F N ■ p(v) = exp(logf N + logp(v)), 
this means that the dependencies of the coefficient of friction as a 
function of logF N will have the same shape for arbitrary velocities, 
only shifted along the logp^-axis by a velocity-dependent shift factor. 
This property gives the possibility to construct dependencies of the 
coefficient of friction on the normal force and the sliding velocity 
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Figure 5 | Measured dependencies of logio// on log w F N at various 
velocities at the temperature 25.5 ± 0.5 °C. 

using a "master curve procedure" similar to those used for 
determining dependencies of the coefficient of friction on velocity 
from measurements at different temperatures 12 : Experimental results 
for the friction coefficient are presented as a function of logF N at 
various velocities in Figure 5. Following this hypothesis, we assume 
that at different velocities, the measured curves are only shifted 
pieces of the same curve. Now, one attempts to shift the curves 
such that they form a single "master curve" (Figure 6). The 
resulting curve gives the dependence of the coefficient of friction 
in a wider range of forces than the range used in the experiment. 
At the same time, the shift factors at different velocities will provide 
the dependence of the coefficient of friction on velocity. The result is 
a complete dependence of the coefficient of friction in a wide range of 
velocities and forces. Repeated for different temperatures and using 
the standard master curve procedure 12 , this will lead to restoring the 
complete law of friction as function of velocity, temperature and 
normal force. However, in the present paper, we avoid the well 
discussed subject of the temperature dependence and concentrate 
our efforts completely on the force dependence. 

Note that the main logic of the result (23) is not dependent on the 
details of the model and even on it dimensionality. The scaling rela- 
tion (23) follows solely from the assumption that the macroscopic 
form of the contact is determined by the macroscopic properties of 
material and do not depend on microscopic details, and on the other 
hand, that the microscopic properties are determined mainly by the 
indentation depth. These general assumptions are equally valid for 
one-, two- and three dimensional models. Below we explain this 
important point in more detail. 

It is well known, that if a rigid body of an arbitrary shape is pressed 
against a homogeneous elastic half-space then the resulting contact 
configuration is only a function of the indentation depth d. At a given 
indentation depth, the contact configuration does not depend on the 
elastic properties of the medium, and it will be the same even for 
indentation of a viscous fluid or of any linearly viscoelastic material. 
This general behavior was recognized by Lee and Radok 22,23 and was 
verified numerically for fractal rough surfaces 15 . Further, the contact 
configuration at a given depth remains approximately invariant for 
media with thin coatings 24 or for multi-layered systems, provided the 
difference of elastic properties of the different layers is not too large 25 . 
It was argued that this is equally valid for media which are hetero- 
geneous in the lateral direction (along the contact plane) 26 . Along 
with the contact configuration, all contact properties including the 
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Figure 6 | Horizontal shifting of the curves shown in Figure 5 relative to 
the curve at the reference velocity of 1 mm/s provides a "master curve". It 
has two distinct linear parts. 

real contact area, the contact length, the contact stiffness, as well as 
the rms value of the surface gradient in the contact area will be 
unambiguous functions of the indentation depth. The indentation 
depth is thus a convenient und robust "governing parameter" for 
contact and frictional properties of media with linear rheology. Note, 
that this is equally valid for tangential contact. This can easily be 
illustrated with the example of contact of a rigid body with an incom- 
pressible elastic half-space: For a circular contact with an arbitrary 
radius a, the ratio of the normal stiffness k z and the tangential stiff- 
ness k x is constant and given by the Cattaneo-Mindlin factor 27,28 , for 
incompressible media kjk x = 1.5. From this follows that for a fric- 
tional contact with the coefficient of friction /.i, the maximum tan- 
gential displacement to the onset of complete sliding is determined 
solely by the indentation depth and is equal to u x> max =1.5 {.id. This 
result does not depend on the form of the body and is valid for 
arbitrary bodies of revolution 29 and even for randomly rough fractal 
surfaces 30 . This fact, that the contact configuration is solely deter- 
mined by the indentation depth is as a matter of fact the only physical 
reason needed to get the simple scaling relations for the coefficient of 
friction between rough rigid bodies and linearly viscoelastic elasto- 
mers described by equation (23). While the particular form (22) can 
depend on the model used, the general functional form (23) is a 
universal one and is not connected with the method of dimension- 
ality reduction used in this paper. 

Experimental. Our numerical and theoretical analysis shows that 
under some conditions the dependencies of the coefficient of 
friction on the normal force, presented in double logarithmic axes, 
is self-similar at different velocities and can be mapped onto each 
other by a simple shifting along the force axis. To prove this 
hypothesis, we measured the coefficient of friction between a band 
of polyurethane (PU) and a steel ball with radius _R = 50 mm. 

The measured coefficients of friction as a function of force are 
shown in Figure 5. If the shifting procedure formulated in last 
Section is valid, all the curves shown in this figure have to be con- 
sidered as different parts of the same curve shifted along the logiv- 
axis. Figure 6 illustrates that it is indeed possible to shift all the curves 
to produce one single "master curve". It is interesting to note that the 
resulting master curve has two distinct linear regions, meaning a 
power dependence of the coefficient of friction on the normal force. 
The crossover between different powers occurs at the force F N ~ 
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Figure 7 | Shifting factors as a function of the sliding velocity. 

40 N which coincides with the force F^=(2/3)G 0 D 3 /«~40N, at 
which the contact diameter 2a becomes equal to the thickness of 
the rubber layer D = 10 -2 m. At this force the stiffness-force depend- 
ence changes from the Hertzian k oc F m to k oc F 1 ' 2 21 , and we expect 
a change of the scaling relation. 

Figure 7 shows the dependence of the shifting factor on the sliding 
velocity. Roughly speaking, the shifting factor is a linear function of 
the logarithm of velocity with the slope —0.15. This result can be 
interpreted as follows. In the intermediate frequency range, the loss 
modulus C'often is a power function of frequency: 

G"(co) = G-(ffl/oj 0 ) /J (24) 

where oj 0 is a reference frequency and fS a power typically in the range 
of 0. 1 to 0.5 see e " 8 " 31 . In this case equation (22) can be rewritten as (here 
for a sphere with radius R, n = 2, c n = 1AR): 



2v/2GVz V 3R 



-i/« 



(25) 



In this case, the shift factor is a linear function of logv with the slope 
— 3/?/2. Comparing this with the experimental value of — 0.15 gives p 
= 0.1. This is compatible both with the rheological data for the used 
rubber compound and with data from literature. 

Discussion 

We analyzed the frictional behavior of elastomers under the follow- 
ing simplifying assumptions: (a) the rigid counter body has a power 
law shape (e.g. paraboloid or cone), (b) the macroscopic contact 
mechanics of the indenter is governed mainly by the low frequency 
shear modulus, which can be assumed to be approximately constant, 
(c) the friction is governed by the corrugations with the smallest 
wavelength in the spectrum of the surface roughness. Under these 
assumptions, we have shown that the coefficient of friction is a func- 
tion of a dimensionless argument, which is a multiplicative function 
of powers of velocity and force. The exact form of this argument 
depends both on the rheology and the macroscopic form of the 
indenter. But independently of the exact form, the dependence of 
the coefficient of friction in the range from very small velocities to the 
plateau occurs to be a universal function of this argument, suggesting 
a generalization of the known "master curve procedure": If the 
dependence of fi on the normal force is presented in double logarith- 
mic coordinates, it will have the same shape for arbitrary velocities, 
only shifted along the velocity axis. We have proven this procedure 
with experimental results obtained on polyurethane rubber. In 



Figure 8 | Experimental set-up for measuring the coefficient of friction. 

combination with the widely used shifting procedure for varying 
temperature 12 , it allows to determine generalized laws of friction as 
functions of velocity, temperature and normal force. The results of 
the present paper generalize and validate the results of the pioneering 
work by Schallamach 32 . 

Methods 

MDR calculation. Method of Dimensionality Reduction (MDR) is based on mapping 
of three-dimensional contact problem to contacts with one -dimensional elastic or 
viscoelstic foundations 18 . It gives exact solutions for contacts of bodies of revolution, 
and provides a good approximation for all properties which depend on the force- 
displacement relationship such as contact stiffness, electrical resistance and thermal 
conductivity, and also dissipated energy and frictional force for elastomers. The 
details are described in the introduction and at the beginning of the Section Results. 

Experiment. The experimental set-up for measuring elastomer friction is shown in 
Figure 8. The rubber band with a size of 300 X 50 X 5 mm was glued to a moving 
stage using a solvent-free two-component epoxy glue. We used polyurethan rubber 
with a static shear modulus of about 3 MPa and yield modulus of 48 MPa. The 
maximum pressure in the contact area was, in all experiments, at least one order of 
magnitude smaller than the latter, so that there was no plastic deformation of rubber. 
The stage could be moved with the aid of a hydraulic actuator with controlled velocity 
in the range of 5 • 10~ 4 m/s to 0.58 m/s. The normal and tangential forces were 
measured with a 3D force sensor, on which the steel ball (radius 25 mm) was 
mounted. The ambient temperature was 25,5' C ( + /— 0,5) and the relative air 
humidity 30% ( + /— 5). Under these conditions the dynamic friction coefficient was 
measured at constant normal force and horizontal velocity. A total of 1680 
measurements were taken. Data for any parameter set (normal force and 
temperature) was averaged over six measurements. Every measurement series was 
started at the smallest normal force, and increased in steps. At every level of normal 
force, the measurement was made with 28 horizontal velocities. Before proceeding to 
the next force level the material was examined for wear both visually and with a 
microscope, and cleaned with pressurized air. At low normal forces wear was virtually 
non-existent, and remained weak even at higher forces. Chemical cleaning agents 
were not used to treat the surface of the rubber. 

Measurements were carried out in the velocity range of 3 ■ 10~ 3 m/s to 2 • 10~ 2 m/s 
for normal forces in the range of 1 N to 1 00 N. The lower velocity bound was chosen 
because the sliding at lower velocities was instationary. The maximum velocity was 
chosen to avoid significant temperature changes in the contact. The local temperature 
rise due to frictional heat can be estimated as AT ~ 2f.tG 0 dv/?., where a ~ 
20 Wm~'K _1 is the thermal conductivity of the steel ball 21 , which almost completely 
controls the thermal flow, and G 0 ~ 3 • 10 6 Pa the static shear modulus of the used 
rubber {Rheological measurements were carried out by Thermoplastics Testing 
Center of UL International TTC GmbH). For the largest force of F N = 10 2 N, we get 

(\ 2 / 3 
3F N /16G 0 R 1/2 J «4-10" 5 m. With fx « 0.5 and v = 

10~ 2 m/s we can estimate the average temperature rise as AT ~ 0.06 T. Maximum 
temperature changes in micro contacts can be estimated as AT ~ 2f.tGoJi 0 v/A, where 
Geo « 1.5 ■ 10 9 Pa is the glass modulus {Rheological measurements were carried out by 
Thermoplastics Testing Center of UL International TTC GmbH) of the used rubber 
and ho ~ 10~ 7 m the rms roughness of the ball, which was determined using a white 
light interferometric microscope. For velocity v = 10" 2 m/s we get an estimation AT 
~ 0.08 K which is of the same order of magnitude as the average temperature rise. 
Due to repeated sliding, the temperature change can get larger than the above 
estimation. The temperature changes of the rubber surface were controlled in 
experiments by an infrared camera (see Figure 8). We found empirically that the 
temperature change does not exceed 1 k for the following range of velocities: up to v 
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~ 4 • 10" 2 m/s for F N ~ 10 2 N, up to v ~ 2 ■ 10" 2 m/s for P N ~ 10 N and up to v ~ 
m/sfor F N ~ 100 N. 
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